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ABSTRACT 

We present a solid relationship between the neutral hydrogen (H I) disk mass and the stellar disk 
mass of late-type galaxies in the local universe. This relationship is derived by comparing the stellar 
disk mass function from the Sloan Digital Sky Survey and the H I mass function from the H I Parkes 
All Sky Survey (HIPASS). We find that the H I mass in late-type galaxies tightly correlates with the 
stellar mass over three orders of magnitude in stellar disk mass. We cross-check our result with that 
obtained from a sample of HIPASS objects for which the stellar mass has been obtained by inner 
kinematics. In addition, we derive the H I versus halo mass relationship and the dependence of all 
the baryonic components in spirals on the host halo mass. These relationships bear the imprint of the 
processes ruling galaxy formation, and highlight the inefficiency of galaxies both in forming stars and 
in retaining their pristine H I gas. 

Subject headings: galaxies: statistics — galaxies: formation — galaxies: evolution 



1. INTRODUCTION 

During the last years ground- and space-based surveys 
allowed to probe the physical properties of many mil- 
lions galaxies both in the local and in the high-redshift 
Universe. These analysis have been mainly focused on in- 
vestigating the stellar component of galaxies, and have 
provided us with a much clearer view of when and where 
star formation occurred along the cosmic time. 

However, there is another baryonic component, namely 
the neutral atomic hydrogen HI, that should be accu- 
rately monitored to understand the process of galaxy 
formation; in fact, such a component just constitutes 
the raw material which stars are made of. According 
to the standard picture, protogalactic halos initially had 
all the same cosmological a mount of HI gas, aro und 1/6 
of the host halo mass (e.g.. iKomatsu et al.ll20lil ). in the 
form of a warm atmosphere. Then a fraction of such 
warm baryons is expected to cool and condense in a cold 
gaseous disk-like component, whereby stars are formed. 
In turn, this cold, star-forming gas can be depleted by 
the energy feedback from type II supernova explosions 
and stellar winds, in an amount modulated by the ratio 
between the total energy injected and the depth of the 
potential well of the host halo; noticeably, the former is 
related to the overall mass of formed stars and hence to 
the galaxy luminosity, while the lat ter crucially depends 
on the mass of the hos t halo (e.g. I White fc ReeslTl978l : 
IFall fc EfstathToul[T980T ). At lower halo masses a large 
HI depletion is due to t he photo-heating by intergalact ic 
UV radiation field (e.g. iHoeft et al.|[2006l: fRicottill2Q09h . 

Therefore, the observational information on the HI 
mass content of galaxies provide crucial constraints on 
galaxy formation theories; a successful scenario must be 
able to reproduce not only the observed stellar mass func- 
tion and luminosity f unction, but also the HI mass func- 
tion (|Mo et alll2005D and the relationships between the 
HI and the stellar/halo mass. 
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Only in recent years, thanks to the completion of rela- 
tively wide blind 21-cm surveys, a wealth of observations 
on HI gas has become available. In detail, Zwaa rTet al.l 
(2005) used the catalog of 4315 extragalactic HI 21-cm 
emission line detections from the HI Parkes All Sky Sur- 
vey (HIPASS) to obtain an accurate measurement of the 
galaxy HI mass function (HIMF) down to an HI mass of 
10 7 - 2 Mq. 

In this work we aim at investigating the relationships 
between the HI mass and two relevant structural prop- 
erties of late-type galaxies: the stellar disk and the 
halo masses. To reach this we exploit: i) a theoretical 
approach that boils down to matching the cumulative 
HIMF mass function and the galactic stellar (or halo) 
mass function; ii) an observational approach that relies 
on a sample of objects for which both the HI and stellar 
disk masses have been directly measured. We show that 
the two approaches agree in indicating a strong correla- 
tion between the gaseous disk and stellar disk (or halo) 
mass. 

The existence, in late-type spirals, of a relationship 
between the HI disk mass and the galaxy luminosity 
is we l l known (e.g. [Robertsl 119751: iRoberts fc Havnesl 
1991 IGavazzi et all Il996h IMcGaueh fc de Blokl Il997t 
Disnev et aiT T2008l so as that between the former 
and the spec tro- photometrically derived mass of the 
stellar disk (IKannapparJ 120041 : IGavazzi et all 120081 : 
iCatinella et all 12010ft . Recently, iCortese et aTT (|201lD 
showed that Hl-to-stellar mass ratio anti correlates with 
stellar mass over ~2 order of magnitudes in stellar mass 
and investigated the effect of the environment on this 
relation. 

Previous results, based on the spectrophotometric es- 
timate of the stellar disk masses, have estabilished the 
existence of the relations subject of the present investi- 
gation, but in a biased way. In fact, especially for spi- 
rals, the luminosity is a poor indicator of the stellar disk 
mass, and, in any c ase, it is uncertain by a factor two 
(jSalucci et al.ll2008h . In addition, it depends on the as- 
sumed initial mass function (IMF) and star formation 
rate, quantities that we would like to study helped by 



2 



Evoli et al. 




12.0 



150 



100 - 



50 - 



5 
Hubble Type 



10 



Fig. 1, — The GSMF of late- type galaxies obta ined from the 
LF by Bcrnardi ct al. (2010) and the M/L ratios by Shankar ct al. 
p00g) is illustrated as a blue line. The GSMF bv lBell et al.1 (12003) 
is also shown (red line) for comparison. 

the results of this paper and not to assume a priori to 
get the results of this paper. Finally, the above relation- 
ships are biased by the fact that spirals with the same 
stellar disk mass, but overabundant or deficient in HI 
content, seem to have different stellar mass-to-light ra- 
tios (and then luminosities) than galaxies "normal" in 
HI content. 

In this work, we aim to estimate the mass of a stellar 
disk, associated to a HI disk in two essentially new, accu- 
rate, model-independent and statistically relevant ways. 
These estimates are expected to yield trustable relation- 
ships or trends, free from biases that are likely to affect 
their interpretation in a cosmological context. Notice 
that Shankar et al.l (|2006f ) by following iSalucci fc Persicl 
(1999), were the first to correlate the kinematical bias- 
free estimates of stellar disk mass with the corresponding 
HI masses, however, their work was based on a sample 
much more limited, in number of objects and magnitude 
extension, than that we use in this work. 

Finally, we apply the cumulative technique to derive 
the relationship between HI and halo masses. Even if 
not strongly motivated as in the previous case, we are 
able to derive a more realistic relationship for these two 
observables with respect to what existing in literature. 

Throughout the paper we adopt the standard value 
Hq = 73 km s _1 Mpc -1 for the Hubble constant, and 
quote uncertainties at 1 — a confidence level. 

2. HI VS. STELLAR MASS RELATIONSHIP 

To investigate the relationship between the stellar and 
the gas mass com ponent in late- type gal axies, we follow 
the procedure by iVale fc Ostrikerl (|2004[ ) and developed 
by I Shankar et al.l ( 20061 ) . First, supported by the evi- 
dence described in section . 1 we assume that, in aver- 
age, the mass of the HI disk is, in statistical sense, an 
(increasing) monotonic function of the mass of the stellar 
disk. 

If two galaxy properties q and p obey a one-to-one 
relationship, we can write: 



<f>(p) dg = ip(q) dq 



(1) 



Fig. 2. — The dis tribution of Hubble - types in the HIPASS galax- 
ies as presented by Rvan-Wcbcr ct al. (2002). 

where ipil) 1S the number density of galaxies with mea- 
sured property between q and q + dq and 4>{p) is the 
corresponding number density for the variable p. The 
solution is based on a numerical scheme imposing that 
the number of galaxies with q above a certain value q 
must be equal to the number of galaxies with p above p, 
i.e., 



(f>{p) dp 



ip(q)dq . 



(2) 



In the following we take p as the HI mass Mhi and 4>(p) 
as the corresponding HIMF, while q as the stellar mass 
M* and tp(l) as the corresponding galactic stellar mass 
function GSMF. T he local HIMF has been measured by 
iZwaan et all (f2 005) u sing th e galaxy data in the HIPASS 
catalog (|Mever et al.l [2004); its shape has been fitted, 
within the range 10 7 - 2 M Q < M m < 10 11 M Q , with a 
Schechter function: 



( f > (M m )dM m = 
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with power law slope a = —1.37 ± 0.03, characteristic 
mass log(MHi/M©) = 9.8 ± 0.03/i^ 5 2 and normalization 
</>hi = (6 ± 0.8) x 10~ 3 ^ 5 Mpc~ 3 dex- 1 . 

Obviously, the two mass functions appearing in Eq. (|2|) 
must be representative of the same galaxy population. 
To check this, we plot in Fig. [2] the Hubble-type dis- 
tribution (obtained from the HyperLeda Catalogue, see 
iPaturel et all [2003| ) of the 1000 bright est HIPASS galax- 
ies as reported in iRvan- Weber et all ([2002. ) . We con- 
clude that the HIMF represents almost entirely disk sys- 
tems: late-type galaxies account for more than 85% (Sb- 
Sc), there is a small contribution from irregular galaxies 
(smaller than 15%), and the contribution from ellipticals 
is negligible (smaller than 2%). 

Thus we calculate the GSMF for late-type and Irreg- 
ular galaxies on the basis of the recent observational re- 
sults reported in IBernardi et al.l (|201CO . Specifically, we 
use their LF for C r < 2.6 (M. Bernardi, private com- 
munication), which implies a small contamination from 

3 C r is the concentration index defined as the ratio of the scale 
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early-type galaxies, around 2% from ellipticals and less 
than 26% from Sa-type objects. From this we build the 
GSMF by adopting the disk Mass-to-Light ratio derived 
from mass modellin g of the (Spiral ) Univ ersal Rotation 
Curve, see Eq. 2 in Shankar ct al. (2006), and we fit it 



with a modified Schechter function (see iBernardi et al.1 
2010, Eq. 9): 



</>(M*)dA/* = <t>* 



T(a//3) Pd \M 



(4) 



with parameters: 0* = 1.05 x 10 -2 Mpc~ 3 , a — 0.385, 
/3 = 0.59 and logM, = 10.05. The function is plot- 
ted in Fig. [1] alongside, for the sake of compa rison, with 
the G SMF of late-type galaxies obtained by iBell et al~l 
(2003) from model-dependent spectrophotometric esti- 
mates of the disk masses. The method we use suffers for 
different uncertainties, in particular in inferring stellar 
masses from kinematical measurements, hence the total 
uncertainties on our results is of the order of 30%. 

More recently the ALFALFA collaboration have pub- 
lished an HIMF based on 10 119 galaxies by pro bing a 
bigger volume than HIPASS (jMartin et aljfeoiOD . The 
new HIMF differs from the HIPASS one at the high mass 
end, to the effect of changing the normalization of the 
Hl-to-stellar mass ratio. We show in Fig. [3] that the dif- 
ferences in assuming the ALFALFA HIMF are within the 
error bars associated to the uncertainties in the HIPASS 
HIMF normalization. 

Then, we solve Eq. ([2]) and derive the relationship be- 
tween the gas to star fraction and the stellar mass; the 
result is shown in Fig. |3] The gas fraction and the stellar 
mass correlates as a broken power-law over about three 
order of magnitudes in stellar mass. Within the mass 
range 10 8 < M* < 10 11 the relationship can be well ap- 
proximated by : 
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(5) 

This relationship, obtained by direct estimate of the 
stellar disk mass, can be compared with that obtained 
by means of the (biased) traditional methods. In Fig. [3] 
we compare our res ult with th e HI t o stellar mass ob- 
tained by using the IBell et al.l (|2003| ) GSMF. The dif- 
ference between the two is particularly pronounced at 
s mall masses , wher e the spectrophotometric M/L ratios 
of lBell et al.l (12003?) are appreciably larger than the kine- 
matical estimates. 

3. HI CONTENT OF INDIVIDUAL GALAXIES 

We derive the relationship between the HI mass (Mm) 
and the stellar disk mass (Md) with a new model- 
independent method by looking at individual late-type 
galaxies. The disk mass is obtained, within a reason- 
able uncertainty, by modeling the galaxy rotation curve, 
whose inner parts are decomposed in halo and disk com- 
ponents. 

Let us first define i? op t = 3.2i?£>, where Rd is the 
exponential thin disc length-scale. This radius, that en- 
closes about 83% of the total light, can be considered the 

which contains 90% of the Petrosian light in the r band, to that 
which contains 50%. 
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Fig. 3. — Ratio of HI to stellar disk mass as function of the latter, 
the dashed area represents the uncertainty related to HIMF nor- 
malization. Triangles represents individual objects (blue symbols 
are for late- types and green for irregulars). Red and orange lines 
show the effects of changing the HIMF and the GSMF as described 
in Sec. 2, respectively. The H2 to stellar mass ratio is also reported 
as a green line. 



physical size of the stellar disc. iPersic fc Saluccil (|1990l ) 
devised a reliable method to estimate the disk mass from 
observational quantities, i.e. from the gravitating mass 
M g inside i? op t (M 9 w G^ 1 V^ pt R op t) and the rotation 
curve logarithmic slope at R op t (V): 



M D = (0.72 - 0.85 V) M g 



(6) 



Wc then proceed to build a sample containing the 
75 objects in HIPASS that have optical photom- 
etry and kinematics of quality sufficient for the 
above method. The r otation curves are t aken from 
I Persic fc Saluccil (119951). lYegorova fc Saluccil (|2007| ) and 
iPri gerio Martins fc Saluccil (j2007l ). By means of Eq.[6]we 
derive the disk mass with a n uncertainty between 10% 
-30% (IPersic fc Saluccilll990l) . 

In Fig. [3] we show the relation obtained for individual 
objects and that obtained by matching the HIMF to the 
GSMF. The two are in very good agreement over two 
order of magnitudes in stellar mass, showing the same 
power law functional form (with slope respectively of 
—0.48 and —0.52) and similar normalization. The agree- 
ment of the individual objects and the statistical relation, 
obtained from two very different methods, indicates that 
the first one is little biased by contamination or incom- 
pleteness of the HIMF, and that the second uses a fair 
sample of individual objects. A stellar disk mass versus 
HI disk mass relation emerges as one of the most impor- 
tant empirical relationships concerning spirals. 

A further gas component in the local galaxies is the 
molecular hydrogen (H2 ) disk. Although we must caveat 
that its mass does not necessarily m onotonically corre- 
late with the stellar disk mass (e.g. ICasoli et al.l 119961 : 
iBoselli et all 120021: iBoker et al.l I2003D . we will proceed 
as above, since more information on this poorly known 
component is certainly needed. 

Let us stress that, unlike the HI mass, the H2 disk 
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Fig. 4. — The baryonic content relative to the initial baryonic 
mass associated to a halo (ft^halo)- Red line refers to stars, green 
line to HI and blue to the total. Dashed line shows the finding by 
IMan'n et aLl 120101') . 

mass estimate relies on indirect tracers as CO lines, with 
uncertai n conversion factors. We adopt the H2MF de- 
rived by lObreschkow fc Rawlingsl ()2009f ) from the local 
CO luminosity function of the CARO Extragalactic CO 
Survey, assuming a variable CO to H2 conversion fac- 
tor fitted to nearby observations. The corresponding 
mass function (H2MF) is well fitted by a Schechter func- 
tion with powerlaw slope a = —1.07, characteristic mass 
log(MH 2 /M ) = 9.2 and normalization </>h 2 = 8.3 x 10~ 3 
Mpc -3 dex -1 . The resulting H2 to stellar mass ratio 
as a function of the stellar mass is shown in Fig. |3j as 
expected, this component turns out to be subdominant 
relative to HI over the whole probed mass range and for 
this reason we do not consider this contribution in the 
rest of the paper. 

4. HI VS. HALO MASS RELATIONSHIP 

It cosmologically relevant to derive the relationship be- 
tween HI mass Mhi and halo mass Afjj in spirals. A 
preliminary step is to obtain the relationship between 
the stellar mass M* and halo mass Mh by the method 
described in Sect. 2. iShankar et al.l (|2006h already ob- 
tained this results but it is worth to redo their analysis 
with updated observational data. To this purpose, we 
need two ingredients. The first is the galactic halo mass 
function (HF ), i.e., the stati s tics o f halos containing one 
single galaxy; IShankar et al.l (|2006f ) evaluated it from the 
standard halo mass function by adding the contribution 
of subhalos, and subtracting the contribution of galaxy 
systems (see their Eq. 9). 

The second is the GSMF of all the local galaxy popu- 
lation, necessary because the HF does not distinguish 
bet ween galaxy m o rphol ogy. We base on the GSMF 
by iBernardi et al.l (|2010l ) without selection criteria in 
concentration index; this is fitted in terms of a modi- 
fied Schechter function (see their Eq. 9) with parameters 
given in their Table B5. 

The relationship derived with these mass functions 
holds for the overall galaxy population, so to proceed 
further we must assume that it also approximately holds 
for each separate Hubble type, in particular, for late-type 
objects. This is justified by the fact that we found that 
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Fig. 5. — In red the HI mass inside the optical radius, in green 
the HI mass outside, and in black the total. 

the fractional amount of the HI component with respect 
to the whole baryonic component vary across Spirals by 3 
orders of magnitudes; on the other hand, from X-ray and 
weak lensing observations, we can infer that galaxies with 
the same halo mass have approximately the same bary- 
onic mass and that, furthermore, the relation between 
the galaxy virial mass and the relativ e baryonic mass 
is ro ughly Hubble Type independent dFukazawa et all 
120061 iNagino fc Matsushitall2009L Ibonato et alJl2009f ). 

Thus, we combine the HI vs. stellar mass relationship 
with the stellar vs. halo mass relationship to obtain the 
HI vs. halo mass relationship. We show the result in 
Fig. [4] the relation can be fitted (to better than 5% com- 
paring with the numerical result) within the mass range 
10 11 M Q < M H < 10 12 - 5 M as: 

Mhi (Mh/9.5 x 10 11 M )°- 33 

9 x 10 9 M Q 1 + (Mh/9.5 x 10 11 M Q )-°- 77 ' 1 ' 

We also plot for compariso n the HI vs. halo m ass rela- 
tionship recently derived by iMarm et al.l (|2010f) by com- 
paring directly the statistics of HI and halo masses. Their 
results appreciably differs from ours since the standard 
halo mass function they adopt includes the contribution 
of galaxy groups systems so it has more objects relative 
to our GHMF; then, the matching procedure of Sect. 2 
leads to a lower HI mass at a given halo mass. 

In Fig. |4] we summarize our results by showing the 
amount of HI and stellar mass (relative to the initial 
baryonic mass )associated with a halo as a function of 
its virial mass. We also plot the overall baryon fraction 
derived by adding the stellar mass to the total gas mass 
obtained by multiplying the HI mass for 1.4 to take into 
account the contribution of He. 

5. DISCUSSION AND CONCLUSIONS 

The correlations of the HI mass with stellar and halo 
masses are extremely relevant in the framework of galaxy 
formation theories. The standard picture envisages that 
every galaxy forms with the same initial amount of 
baryon in the form of HI gas, and what we observe now 
is the left-over of the processes that took place during 
galaxy formation. 
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Fig. 0] shows that late- type galaxies are extremely inef- 
ficient in retaining their initial baryon content, i.e., most 
of the initial HI gas has been removed from the host halo. 
Less than 10% is retained in galaxies with halo masses 
below 10 11 M Q , and this value drops to few % for halo 
masses above 10 12 M@. Such a behavior is likely due to 
supernova feedback. Thus only a small fraction of the 
initial baryon content is eventually exploited for the star 
formation. Note that, in massive halos, stars are the 
dominant baryonic component while in smaller halos HI 
gas is. 

Let us stress however that the baryon cycle in spirals is 
very complicated to understand. It may depend, in addi- 
tion to SN feedback, IGM ionization, gas cooling time, on 
the interplay between galaxies and their environment, es- 
pecially for low mass halos. Notice SPH simu lations have 
not ye t converged to a definitive result, e.g. IHoeft et al.l 
(|2006() finds that halos with M > 10 10 5 M fl a r e able 
to retain all their baryons while [Pilkington et al.l (|201 lh 
find the galaxy formation process able to remove most of 
the original baryonic material. All this means that the 
processes that are responsible of the evolution of galactic 
gas about which this paper provides valuable information 
are not fully understood. 

One can wonder why this gas, although not being 
ejected by supernova feedback, has not been used for 
star formation. To answer the question, we look at where 
this residual HI gas is presently residing by highlighting 
in the previous correlations the contribution from the HI 
gas located inside or outside the stellar disc radius i? op t- 
Therefore, we model the gas surface density of late-type 
galaxies wit h the functional fo rm observed in most Spi- 
rals (see e.g. iBigiel et al.H201Q[) 

10g ^ " \ log So - 2 (r - R op t)/(Ri - Ro P t) if r > R op t 

where R{ is the radius at which the surface density drops 
at 1/100 of the value at R op t that we assume as the size 
of the HI disk and So is the HI surface density central 
value. 

We need now to relate the lenght-scale of the stellar 
distribution with that of the neutral gas. Notice that 
our aim is to obtain qualitative results, in this view the 
assum ptio ns we take are well justified . IBroeils &: Rhed 
(|1997D and lRhee fc van Albadal (| 19961) published the HI 
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surface density profiles for 60 spirals of known optical 
radii i? op t and blue luminosity Lb ( that are given in 
Tab. 1 of iRhee fc van Albadal (jl996[ ). notice that the 
quantity in the fourth column is ~ i? pt/2). From these 
measurements they derived: a) the HI half-mass radius 
R e g', b) the total mass Mhi (given in columns 5 and 6 
of the same Table). From these quantities we obtain 
a strong R e ff vs i? op t relationship and, by the defini- 
tion of R e s, the relationship Rf — F(R c g(R opt ), R op t)- 
Moreover, to transform light in stellar mass, we use, 
without loss of gene rality, M* = 10 n (L B /l0 11 L Q ) 1A M Q 
(jSalucci et al.l f2007). Finally by combining and manipu- 
lating the above empirical relationships (that also imply 
to assume Eq. 8) we obtain: 

Rf/Ro P t = 3-2/3 log (M*/10 9 M Q ) . (9) 

The above indicate, not surpisingly , that small galaxies 
have a larger HI disk, in terms of the stellar disk size. In 
Fig. [5] we show how the HI mass is divided in those inside 
and outside i? pt, the radius inside which the stars reside. 
The former is the dominant component for massive ob- 
jects, while the latter gives a dominant contribution in 
small galaxies. 

The overabundance of HI over stellar mass in small 
objects, is due to material located far away the stellar 
disc and mostly unprocessed. It is worth noticing that in 
these objects at these radii the HI surface density is much 
lower than the threshold of order 1 Mq kpc~ 2 needed by 
the Toomre criterion to form stars. This HI component 
has not been at disposal for the latter process and it 
never will. Let us stress that the inefficiency of the star 
formati on process in the o uter regions of discs is directly 
probed IBigiel et al.1 (pOlOl ). 

To sum up, in this work we have derived robust correla- 
tions between the HI and stellar (halo) mass for late-type 
galaxies in the local Universe. These relationships bear 
the imprint of th e processes ruling galaxy formation (see 
ICook et al.ll2010l for a theoretical approach that consider 
them), and highlight the inefficiency of galaxies both in 
forming stars and in retaining their pristine HI gas. 
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